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Abstract 

Our Galaxy is the largest nuclear interaction experiment which we know, because of 



Q^ ■ the interaction between cosmic ray particles and the interstellar material. Cosmic rays are 

o : 

^ ■ particles, which have been accelerated in the Galaxy or in extragalactic space. Cosmic rays 

^ ■ come as protons, electrons, heavier nuclei, and their antiparticles. Up to energies up to 

O ' some tens of TeV of particle energy it is possible to derive chemical abundances of cosmic 

^ ■ rays. It has been proposed that cosmic ray particles can be attributed to three main sites 

of origin and acceleration, a) supernova shocks in the interstellar medium, b) supernova 
;_j ■ shocks in a stellar wind of the predecessor star, and c) powerful radio galaxies. This 

proposal leads to quantitative tests, which are encouraging sofar. Quantitative models for 
transport and interaction appear to be consistent with the data. Li, Be, B are secondary 
in cosmic rays, as are many of the odd-Z elements, as well as the sub-Fe elements. At 
very low energies, cosmic ray particles are subject to ionization losses, which produce a 
steep low energy cutoff; all particles below the cutoff are moved into the thermal material 
population, and the particles above it remain as cosmic rays. This then changes the 
chemical abundances in the interstellar medium, and is a dominant process for many 
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isotopes of Li, Be, B. With a quantitative theory for the origin of cosmic rays proposed, 
it appears worthwhile to search for yet better spallation cross sections, especially near 
threshold. With such an improved set of cross sections, the theory of the interstellar 
medium and its chemical abundances, both in thermal and in energetic particles, could 
be taken a large step forward. 



1 Introduction 



Before we can use our Galaxy as a tool for the interaction of cosmic rays and thermal 
material, we need to understand the origin of cosmic rays. The origin of cosmic rays is still 
a question [^l], ^ P] which is not finally answered; however, already some 



time ago Cocconi pSf argued convincingly that the very high energy cosmic rays must 
originate outside our Galactic disk, since their Larmor motion could not be contained. 
While the questions about the subtleties of cosmic ray acceleration provide ample material 
for discussion, the debate about the origin of cosmic rays of moderate energy has reached 



a consensus, that they are produced in the Shockwaves of supernova explosions 0, 
|48| , |2^ , |T6|, ^, be it into the interstellar medium, or into a stellar wind 

j, However, the origin of the cosmic rays of the highest energy has remained 
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under dispute. Many of the relevant issues here have been dealt with in the excellent 
review by Hillas (1984 [Q) and in the books by Berezinsky et al. (1990, 0) and Gaisser 
(1990, 0). 

Here we are concerned with the interactions of cosmic rays in the Galaxy, and so 
we will adopt the picture that indeed the cosmic ray particles originate in the shocks of 
supernova explosions. 
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Using this concept (see, e.g. the review by Ginzburg |^), we will describe recent 
advances in our theoretical attempt to formulate a quantitative theory for the cosmic rays 
in the Galaxy. The interaction between energetic particles and the interstellar medium 



has long been of interest p8| , p^ . We observe consequences of such interaction, such as 
gamma ray emission in lines or in continuum, as well as abundances of some elements and 
isotopes (see the comprehensive review by Reeves and the account given by Bloemen 



[|18| ). A recent example of a new measurement of the Boron isotope ratio, together with 



a summary of relevant references, has been given in The detection of gamma ray 

lines, presumably from excited nuclei after nuclear collisions between energetic particles 
and interstellar medium nuclei (predicted a long time ago by Meneguzzi & Reeves ||50|| , 



and Ramaty et al. [^), from the Orion complex [0] has aroused the interest of many 



[p!9| , [20| , |54| , [79[| . Especially the group around R. Ramaty has contributed to the 



discussion, based on their experience with energetic particle interactions in the solar 



activity regions [6^, ^ |6^, 0. The situation has possibly improved, as 

we will try to demonstrate, since we have now a quantitative proposal to account for the 
origin of cosmic rays, and while many of the aspects of this proposal remain to be worked 
out and verified, it may provide a useful basis for further investigations. Therefore here 
we will try to demonstrate that it will be worthwile to obtain better cross sections for 
many of these interactions, so that these interactions may become a quantitative tool in 
the future. 

The structure of this review is as follows: First we briefly summarize the recent pro- 
posal to account for the origin of cosmic rays; then we describe some aspects of injection 
of cosmic rays, and their electromagnetic interaction with the interstellar medium gas; 
then we go through the arguments for the various interaction sites, near the source and 



far from the source; for the latter argument we go through the concept of trapping and 
leakage from interstellar clouds in some detail, since it is new. Finally we draw some 
conclusions and stress the importance of better cross sections. 



2 A quantitative proposal for the origin of galactic 
cosmic rays 

Cosmic rays arrive at earth with energies from several hundred MeV/particle to 310^'^ 
eV; their spectrum for protons is at GeV energies close to i?"^ ''^, and for He and higher 
elements close to i5~2.65 i^gjQ-^ ^ knee at ^ 5 10^^ eV, where the spectrum turns down 
to about E^^'^, to flatten out again near 3 10^*^ eV, called the ankle {e.g. p9| , |89[|). 
The chemical composition is roughly similar to that of the interstellar medium, with 
reduced hydrogen and helium relative to silicon, and with the same general enhancement 
of elements of low first ionization potential as we find in solar energetic particles. The 
low energy end of the observed spectrum is cut off due to interaction with the solar wind. 
There is reason to believe that in interstellar space the cosmic ray spectrum extends far 
below what we can observe at Earth. 

In the newly proposed theory (starting with |^) the origin of the cosmic rays below 

3 10^^ eV is traced to the Shockwaves caused by supernovae exploding either into the 
interstellar medium, or into the predecessor stellar wind, following some rather classical 
ideas; the new element is a premise on the particle transport in the shock region, inspired 
by the observations of the radio polarization in supernova remnants, and the actual motion 
of radio features, as well as the size of the observed X-ray and radio supernova remnant 



shells [13|: These data suggest a strongly turbulent interaction region rather than a 



smooth shock wave, consistent with several arguments which have demonstrated that 



cosmic ray influenced shocks are unstable (see ^ and the detailed discussion of this 



point in [|T^]). This premise is the principle of the smallest dominant scale, which follows 



work by Prandtl (1925 |5^) and von Karman & Howarth (1938 This principle is 

used to find a length scale and a velocity scale, describing turbulent transport. Applied to 
supernova shock shells, this principle leads to some fraction of the radius of the spherical 
shock as a length scale and the velocity difference across the shock as the velocity scale, 
associated with fast convective shock turbulence, and therefore to a specific model of the 
transport of particles in the shock region. In the construction of a transport coefficient for 
energetic particles, then these scales are used, and thus determine, e.g., the time which a 
particle spends on either side of a shock; this time scale is in turn important for adiabatic 
losses, which a particle experiences, as well as energy gains by drifts in the electric fields, 
seen in the moving shock frame, and thus determines the spectrum of the final particle 
spectrum. This then gives net an appreciable energy loss during the expansion of the 
supernova shock, and leads to a steepening of the predicted spectrum as compared to the 
plane-parallel shock case. 

Figure 1. A schematic picture of the proposed three different source sites and their 
respective contributions (adapted from |J7^). There is a contribution from supernovae 



exploding into the interstellar medium, component 1. The next two components arise 
from supernovae exploding into a predecessor stellar wind, components 2 and 3; the polar 
cap contribution, 3, comes from the polar region of the acceleration in wind-supernovae. 
Finally, component 4 comes from the hot spots of radio galaxies. 



The proposal leads to quantitative predictions for i) the spectra both below and above 
the knee of the cosmic ray spectrum near 5 10^^ eV, where the spectrum turns downwards, 
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ii) the particle energies of the knee and the various cutoffs, as well as iii) the chemical 
composition. We have been able to subject these predictions 0, ||, Q to a variety of 



tests in various publications {e.g. [|78|) and reviewed them as well; the latest overviews of 
these developments are |Ty, |13|, |T^. We continue to perform further tests using ever more 
detailed and newer data. 

2.1 Summary of the predictions for nuclei 

The proposal is that three sites of origin account for the cosmic rays observed, i) supernova 
explosions into the interstellar medium, ISM-SN, ii) supernova explosions into the stellar 
wind of the predecessor star, wind-SN, and iii) radio galaxy hot spots. Here the cosmic 
rays attributed to supernova-shocks in stellar winds, wind-SN, produce an important 
contribution at all energies up to 3 10^ GeV. 

Particle energies go up to 100 Z TeV for ISM-SN, and to 100 Z PeV with a bend at 600 
Z TeV for wind-SN. Radiogalaxy hot spots contribute up to about 100 EeV at the source, 
with some sources up to 4 ZeV, = 410^^ eV [Q. These numerical values are estimates 
with uncertainties of surely larger than a factor of 2, since they derive from an estimated 
strength of the magnetic field, and estimated values of the effective shock velocity. 

The spectra are predicted to be E-2-75±o.o4 ISM-SN, and ^-2.67-o.o2±o.o2 
SN below the knee, and i^^3.07-o.o7±o.o7 ^^^^ wind-SM above the knee, and E""^'^ at injection 
for radiogalaxy hot spots. The polar cap of the wind-SN contributes an £'~2-33 component 
(allowing for leakage from the Galaxy), which, however, contributes significantly only near 
and below the knee, if at all. These spectra are for nuclei and are corrected for leakage 
from the galaxy. 

The chemical abundances are near normal for the injection from ISM-SN, and are 



6 



strongly enriched for the contributions from wind-SN. At the knee the spectrum bends 
downwards at a given rigidity, and so the heavier elements bend downwards at higher 
energy per particle. Thus beyond the knee the medium nuclear mass elements dominate 
all the way to the switchover to the extragalactic component, which is, once again, mostly 
hydrogen and helium, corresponding to what is expected to be contributed from the 
interstellar medium of a radiogalaxy, as well as from any intergalactic contribution mixed 
in P]. This continuous mix in the chemical composition at the knee already renders 
the overall knee feature in a spectrum in energy per particle unavoidably quite smooth, 
a tendency which can only partially be offset by the possible polar cap contribution, 
since that component also is strongest at the same rigidity, where the bend in the overall 
spectrum occurs; this term rigidity refers to the factor occurring in the expression for 
the Larmor radius for any energetic particle, and stands for p/Z, the momentum divided 
by the charge; thus nuclei at the same rigidity have the same Larmor radius in their 
gyromotion in a magnetic field. 

2.2 Observational tests 

These predictions can be compared at some detail with data, and we have given adequate 
comparisons in previous work; a summary of the predictions and tests is given in Table 



1, adapted from |11]: 



Table 1. Spectral indices for hydrogen, helium and heavier nuclei. 
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Experiment 



Energy Range 



element range 



sp. index 



Predicted 



Webber 
LEAP J72 
JACEE 
JACEE 
Sokol 



m 



Ryan et al. 
MSU ^1 



JACEE i 
Japan 

predicted 



LEAP J72 
RICH 



Ryan et al. |71| 
Sokol 



JACEE i, I 
Japan [^| 
Sokol II 
Akeno |5l]| 
Akeno ||7^ 
Tibet AS7 i 

predicted 



HP [|9 



HP P 



FE P] 
Akeno 
Akeno |51[ 
Tibet AS7 

FE ra 



Akeno fSTJ 



below knee 
1-50 GeV 
10-100 GeV 
<40 TeV 
>40 TeV 
>5 TeV 
50-2000 GeV 
10-200 TeV 
50-200 TeV 
8-50 TeV 



below knee 
10-100 GeV 
100-1000 GV 
50-2000 GeV 
>5 TeV 
50-200 TeV 
8-50 TeV 
>5 TeV 

< 5 10^5 eV 
below knee 

< 101^-^5 eV 



above knee 

< 0.410^^ eV 

0.4 - 4 lO^s eV 

2 - 4 10^^ eV 

above knee 

5 10^^ eV - 6 10^^ eV 

> W^-^^ eV 

2 10^^ - 4 10^9 eV 

610^^-710^^ eV 



H 2.75 ±0.04 

H + He 2.70 ±0.05 

H 2.74 ±0.02 

H 2.64 ±0.12 

H 3.22 ±0.28 

H 2.85 ±0.14 

H 2.75 ±0.03 

H 3.14 ±0.08 

H 2.77 ±0.06 

H 2.82 ±0.13 



He,..,Fe 
He 
He 
He 
He 
He 
He 
all 
all 
all 
all 



all 
all 
all 
all 
all 
all 
all 
all 
all 



2.67 ±0.02 

2.68 ±0.03 
2.64 ±0.09 
2.77 ±0.05 
2.64 ±0.12 

2.67 ±0.08 
2.75 ±0.15 

2.68 ±0.07 
2.62 ±0.12 
2.66 ±syst. 
2.60 ±0.04 



3.07 ±0.07 

3.01 ±0.02 
3.14 ±0.06 
3.07 ±0.01 
3.07 ±syst. 

3.02 ±0.03 
3.00 ±0.05 
3.18 ±0.01 
3.18 ±0.08 



±0.02 



±0.07 



We note that the error distribution of the prediction below and above the knee, for the 



elements He and higher nuclei is asymmetric with respect to the central prediction. The 



systematic errors inherent in the analysis given in ^7j, and indicated as such in the Table, 
cannot be easily quantified, since they arise from the errors in the Monte-Carlo used for 
the modelling the airshowers; however, the fit to the data is quite acceptable, and so we 
believe that this systematic error is small. The cutoffs in the three source components and 
their chemical abundances can be checked using vertical and slanted airshowers, and are 
all consistent to within 20 % with prediction [^. The gradual switch from one spectral 
range to another across the knee is clearly recognizable for the Tibet AS7-experiment, 
for which this energy range is about a factor of 10, consistent with the expected gradual 



change in chemical composition (see [T^l)- The last two lines in the Table refer 



to energy ranges which cover some of the ankle, where the spectrum varies, due to the 
switch to a new contributor, the expected extragalactic cosmic rays. Here we note also 
that the cosmic ray spectra of the various chemical elements and electrons can be studied 



separately, and all are consistent with the predictions in the GeV to TeV range ^ 
This is the range of interest here. 

3 Injection of cosmic ray nuclei 

For the elements He,..C, O,.. Fe the injection law can be written as a powerlaw in 
momentum p 



N{p) ~ p-'-'' dp, (1) 

which extends all the way down to non-relativistic energies. This means that with p = 
AuipCy (3, where A is the atomic weight of the nucleus considered, and 7 and /3 are 
the Lorentz-factor and velocity in units of the velocity of light c, the spectrum at sub- 
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relativistic energies can be written as ~ 13'"^'^^ d(3. 

The energy loss in interactions with electrons, unbound (then proportional to ne, the 
density of free electrons) or bound in a shell around a hydrogen nucleus (then proportional 
to the density of neutral hydrogen atoms; heavier elements can normally be neglected 
here) of the thermal matter can be written as 

where Z is the charge of the energetic nucleus losing energy. This simple behaviour is 
valid only for suprathermal energies and sub-relativistic speeds (see, e.g. ||52|| ). 

After traversal of thermal matter for some time r the interaction results in a low 
energy cutoff of the distribution of energetic nuclei, and a law of 

~ [5^d[5 (3) 
below the cutoff, and the original law above the cutoff. The cutoff energy is given by 

~ {Z\ne,nH)TY'^. (4) 

All the particles which are lost to the energetic particle spectrum below the cutoff are 
shifted in phase space to the thermal particles, and can modify the chemical abundances 
there. This effect is especially important in the case that the chemical abundances in 
energetic particles are very different from those in the interstellar medium, and this is the 
case for some elements, such as for many isotopes of Li, Be, B. 

The column density along the twisted and scattering path of a charged particle in a 
highly chaotic magnetic field configuration is referred to as grammage, and this grammage 
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is the relevant quantitity to discuss cosmic ray interactions. This grammage can be 
inserted into the above expression, and then leads to estimates of the cutoff energies of 
order 100 MeV for hydrogen and correspondingly more for heavier nuclei. 



4 Spallation of cosmic ray nuclei 

Cosmic ray nuclei can be broken up in collisions with thermal matter; this process is called 
spallation. Obviously, there is a corresponding interaction between energetic protons, 
and thermal material comprising heavier nuclei such as Carbon. In such collisions the 
remaining nuclei can also be excited, and then emit 7-ray lines. 

There are several sites, which can be distinguished, where spallation is relevant (see. 



e.g.^ the recent work in this area |3^, Q ^ ^ 
4.1 Sites of spallation 

First of all, the massive stars, which explode as supernovae after going through an evo- 
lutionary phase accompanied by heavy mass loss, usually have a molecular cloud shell 
around their wind-zone. When the central star explodes, it gives rise to a powerful shock 
wave, which races through the wind, and then smashes into the shell [Q; since the shock 



is loaded with energetic particles, these particles then spallate in the shell. From the 
abundance of sub-Fe elements one can estimate that the grammage in this shell is of 
order Ig/cm^ |8^, consistent with the data from radio and millimeter observations. 
This apparently is the dominant process at higher energy to account for the abundances 
in cosmic rays for most odd-Z elements, for the sub-Fe elements, and for some Li, Be, and 
B isotopes. 

In this case the spectrum of the secondary particles A^^ is the same as the primary 
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particles Np. 



Ns ~ Np. (5) 

Next is the interaction in clouds, and here we have to distinguish between the energy 
range for which the particles move diffusively through a cloud, and the higher energy 
range, where they move unencumbered through the cloud material. It is this latter ap- 
proximation which is commonly used in the literature. 

The secondary particles are then created in the clouds, and diffuse out of the galaxy, 
and so their creation equation can be written as 

dNs N„ Ns 

-fT = — (6) 

at Ts TL,gal 

where is the spallation time scale, and r^^gai is the time scale for diffusion out from 
the disk of the Galaxy. There is a fair amount of evidence that this latter diffusive 
transport can be derived from a Kolmogorov spectrum of interstellar turbulence []70| , p9| . 
The evidence for such a law of turbulence in the ISM has been discussed extensively in 



12, B3]. The solution to this equation is in the stationary case 



Ns = Np^, (7) 

which translates to an energy dependence of the ratio of secondary to primary isotopes 
and elements of 
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in the case of a Kolmogorov spectrum; here we have neglected for didactic simphcity the 
energy dependence of the spallation. The ratio of secondary to primary nuclei has been 
used in the past to argue that in fact the spectrum of interstellar turbulence is not a 
Kolmogorov law. Since the boron/carbon ratio B/C gives an energy dependence of close 
to E^^-^ a Kolmogorov law did not seem to be consistent with the data. 

However, this line of argument is only true, if the cloud interaction is stationary; on 
the other hand we do know that interstellar clouds have their own temporal evolution, 
and so we need to check what happens when clouds form and dissipate again, e.g. by 
heating from newly formed stars. The decisive difference to the argument above arises, 
when we consider the formation of clouds, and we will proceed to do this in the next 
section. 

4.2 The capture of cosmic rays in clouds 

Here we wish to explore the following concept: The interstellar medium is forming large 
molecular clouds out of its small fragments and warmer parts. Gravitational instability is 
a key process in the collapse of clouds or cloud fragments. Gravitational instability sets 
in, as soon as the time scale for free-free collapse is shorter than the time scale for any 
pressure signal to propagate through the cloud. This means that the collapse also needs 
to be faster than the Alfven velocity [|76[. As a consequence, cosmic rays are trapped 
upon the formation of a gravitationally bound system, such as a molecular cloud, since 
cosmic rays cannot stream significantly faster than the Alfven velocity [Q. 

Trapped cosmic rays can get out of the cloud by diffusion; diffusion is a good approx- 
imation only as long as the mean free path for scattering by magnetic irregularities is 
significantly shorter than the size of the cloud. This entails an upper particle energy limit 
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for the diffusion approximation. 

Consider tfien a particle population of cosmic rays Np^i{E, t) trapped in a cloud, where 
the index 1 stands for inside: 



dt TLrl 



with 



TL,cl = TlM^)-'''- (10) 
-C/Q 

This energy dependence follows from the concept that small scale turbulence in media, 
which are magnetic and partially or fully ionized, can be approximated by a Kolmogorov 

law, as discussed above. 
The solution is clearly 



Np,i = Np,i,o{E)exp{—^). (11) 

'Tl,cI 

The particle population outside the cloud, but coming from inside, is then given by 



dNp,2 , iVp,i 



a t tl,cI 

which translates to 

Np,2 = Np,,,,{E){l - exp{-^^^l^^)} (13) 

TL,clfi 

Secondaries are produced in nucleus-nucleus collisions inside the cloud, and so their 
production equation reads 
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dt Ts Tl,cI 

The solution is 

N,,r{E) = Np^^,^{E)-exp{—^). (15) 

The secondaries outside the cloud are just those produced inside and leaking out, and 
so we have the relation 



= H — ^. (16) 
at tl,ci 

The solution to this differential equation is then 

NsAE) = ^^^^r^,d fx'e-^'dx', (17) 

Ts Jo 

where x = t/rL^ci- This entails for long times then 

NsAE) = ^^^r,,iA^r^'- (18) 
Therefore, the secondary particles, injected into the interstellar medium outside the 

original cloud, have a spectrum which is steeper than the primary particles by 1/3. Or, 

given that the primary particles are well approximated by a spectrum of E~^/^, the 

secondary particles at injection have a spectrum of E~^. 

Now, considering then also the leakage from the Galaxy generahzes this results and 

gives the equilibrium spectrum for secondaries: 



dNs^ ^ _ ^^g^ 

dt Tl,cI rL,gal 
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The solution then is 



= TL,ci exp{ —) (1 - rL,ci/rL,gair'' r x'e-^'dx'. (20) 

Ts TL,gal JO 

where now 

x = t{-^--^). (21) 

iL,cl 'L,gal 

Without loss of generality we can assume that tl,cI < TL,gai, when the integral con- 
verges; in the opposite case a brief calculation confirms also the convergence. 

The next step is to assume that we are at present at no particular time; for each 
individual source this corresponds of Ns^2 to an integration over past injection time to 
give A^*2j ^^^^ many sources then no longer changes the spectrum, but only the 

normalization. Clearly, after some long time, the remnants of the cloud are dispersed, 
but then the residual population of secondaries is no longer significant; this then ensures 
that the sum over many sources does not diverge. This then means the further integral 
gives already the proper energy dependence 

= ^^^^ ^lAE) tlUE) (1 - r,,ci/rL,gai)-' m, (22) 

with 

I(t) = r e"^' dx' r se-'ds, (23) 
Jo Jo 

where 

X'' = t/TL,gal, 
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^ = x'TL,gali^/rL,d-^/rL,gal)- (24) 

The ratio of x/x' is energy independent, since in our concept both the diffusion from 
the cloud and the leakage from the Galaxy have the same energy dependence. The integral 
I{t) can be worked out in closed form analytically, and approaches a constant value for 
reasonable large times x" ^ 1. 

The energy dependence of the secondaries, as compared to the primaries is then clearly 

N,,2/Np,i ~ E-'/\ (25) 

with our modelling of the interstellar and intracloud turbulence with a Kolmogorov spec- 
trum. 

This is in accord with the observations, such as by Engelmann et al. This is 

in contrast to the usual finding that a stationary leaky box gives a ratio of secondary to 
primaries ~ if we use a Kolmogorov spectrum for turbulence. 

Therefore, considering the non-stationarity of the normal interstellar medium, we can 
readily explain the ratio of secondaries to primaries, and at the same time use a spectrum 
of turbulence which is consistent with all other observational evidence. Key was the use 
of the gravitational instability condition for the formation of a cloud. 

Translating this result into the language common in the literature, this means that 
escape length as measured in gm/cm^ and escape time can no longer used synonymously. 
The escape time is given by r^^gai, and is proportional to E~^^^ in the relativistic range of 
particle energies. The escape length as a means to describe interaction has three different 
regimes, and the one relevant in the GeV/nucleon range is, as before, about E^^'^, and 
here, in our simplistic model, ~ E^"^^'^. 
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4.3 Energy dependence of secondary /primary ratio 

In the following we adopt for nuclei such as He and higher in mass the primary cosmic 
ray spectrum of E~'^^^ at injection and the Kolmogorov law of turbulence, giving an 
energy dependence of a diffusive time of E~^/^. Therefore, the energy dependence of the 
secondary to primary ratio has three simple domains, which can be summarized as follows: 

• The spallation in the molecular cloud shell around exploding massive star winds 
leads to a ratio of secondary to primary nuclei as a function of energy in the inter- 
stellar medium observable of 

Ns/Np ~ canst. (26) 

• The spallation in the energy range where trapping occurs for cosmic ray nuclei leads 
to 

Ns/Np ~ E-"^'^. (27) 

• And the higher energy range when the interaction is no longer diffusive, we return 

to the canonical solution, which in our case gives 

Ns/Np - E~^/^. (28) 

A comparison with the data suggests that we discern only regime 1 and 2, and that 
regime 3 is never a dominant contributor. The data suggest that the switch between 
regime 1 and 2 occurs near an energy per nucleon of about 20 GeV/n. To repeat, the 
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spallation is described by the first two domains given above, and the escape time corre- 
sponds to a E~^/^ law throughout the relativistic particle energy range. 

5 Chemical abundances 

The origin of the chemical elements and their isotopes can be traced to three main source 
sites (see K 



The big bang nucleosynthesis accounts readily for H, ^He, ^H, ^He, and ^Li. Deu- 
terium, after some excitement about absorption lines in quasars, seems to be now 



in agreement given the first measurements in a neighboring galaxy ||22|. Thus big 
bang nucleosynthesis does seem to give a coherent picture of a universe, where only 
a small fraction of the critical density is made up of normal baryons. 

Stellar interiors and stellar envelopes provide clearly most heavy elements, spewed 
into interstellar space in supernova explosions. Some light isotopes such as deu- 
terium are destroyed in the interior of stars. 

The interactions of cosmic rays with thermal matter can explain a number of features 
both in the abundance distribution of thermal matter, as well as in the distribution 
of cosmic rays: First, the even-odd-Z distribution is dissimilar between the interstel- 
lar medium and the higher energy cosmic rays, with spallation providing a higher 
abundance for the odd-Z elements of cosmic rays. Second, the sub-Fe elements in 
the cosmic rays are also due to spallation. And finally, most isotopes of the elements 
Li, Be, and B are provided by cosmic ray interaction both in the interstellar medium 
and in the cosmic rays. 
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One test |^2|, ^] is the effect of ionization losses on the low energy protons, which 
provide also an ionization and heating source in molecular clouds; it is an important test 
for the entire concept that the cutoff in the proton spectrum due to such losses is consistent 
with the cutoff in the spallation product spectrum required to explain the abundances of 
Li, Be, and B in the interstellar medium. This is the case. 

There is a large amount of work yet to be done, to test the detailed concept proposed, 
in order to account for the chemical abundances in some detail, for the abundances of 
radioactive isotopes, and for accurate isotope ratios. This will provide stringent tests for 
this theory as for any other, and may yet disprove it. 

6 Outlook 

Given that a quantitative theory is beginning to show the promise of an explanation for 
the origin of cosmic rays, it may be worthwile to obtain much better cross sections for the 
cosmic ray interactions, especially near the critical threshold for any reaction. This would 
then allow to not only provide a quantitative explanation of the various abundances, but 
also to actually use them to study both cosmic rays and the interstellar medium. 
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man Nath, Ray Protheroe, Reuven Ramaty, Wolfgang Rhode, Marcia and George Rieke, 
Dongsu Ryu, Eun-Suk Seo, Todor Stanev, Alan Watson, and Barbara Wiebel-Sooth. The 
new element, the concept of time-dependent trapping in interstellar clouds, was devel- 
oped during the Calgary conference 1993, and then further evolved in many discussions; 
I wish to thank the organizers of the Calgary conference and also the organizers of many 
subsequent conferences for generously inviting me; the meeting MESON96 at Krakow 
1996 - where this lecture was given - has been especially stimulating. I thank all my 
discussion partners and apologize for any errors and omissions which surely remain in the 
manuscript. 
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